Magnetic double gradient instability and flapping waves in a current sheet 
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A new kind of magnetohydrodynamic instability and waves are analyzed for a current sheet in 
the presence of a small normal magnetic field component varying along the sheet. These waves and 
instability are related to existence of two gradients of the tangential (B-r) and normal (-B„) magnetic 
field components along the normal (VnBr) and tangential (Vr-Bn) directions with respect to the 
current sheet. The current sheet can be stable or unstable if the multiplication of two magnetic 
gradients is positive or negative. In the stable region, the "kink"-like wave mode is interpreted as 
so called flapping waves observed in the Earth's magnetotail current sheet. The "kink" wave group 
velocity estimated for the Earth's current sheet is of the order of a few tens kilometers per second. 
This is in good agreement with the observations of the flapping motions of the magnetotail current 
sheet. 
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INTRODUCTION 

Thin current layers are typical structures in the Helio- 
sphere, including the solar corona, solar wind and plane- 
tary magnetospheres. We address some of the magneto- 
hydrodynamic aspects concerning the stability of current 
layers which are still poorly understood. In particular, 
CLUSTER observations in the Earth's magnetotail cur- 
rent sheet indicated the appearance of strong wave per- 
turbations propagating across the current sheet. Many 
event studies indicated very large current sheet varia- 
tions and a predominant wave propagation in the trans- 
verse direction with respect to the magnetic field plane. 
The existence of such kind of waves associated with flap- 
ping motions was confirmed in many statistical studies 

E B, B i, i, IS 01 which allowed one to identify them 

as the "kink"-like perturbations. The plasma sheet flap- 
ping observations are interpreted as crossings of a quasi- 
periodic dynamical structure produced by almost vertical 
slippage motion of the neighboring magnetic flux tubes. 
The frequency of the flapping motions, estimated from 
observations is Wf'^ 0.035 s~^ [if. For a majority of the 
observed events [J| , a group speed of the flapping waves 
was found to be in the range of a few tens (30-70) kilome- 
ters per second. The wavelengths and spatial amplitudes 
are estimated to be of the order of 2 - 5 Re {Re is the 
Earth's radius) 0] • 

A preferential appearance of one ("kink" -like) mode of 
the flapping motion was reported by 0]. CLUSTER ob- 
servations give rise to the assumption that the flapping 
motions are notably more frequent in the central part of 
the tail than near the flanks. In the near-flank tail regions 
the motions of flapping waves are predominantly from 
the center to the flanks [2|. These experimental results 



conflrm an internal origin of the flapping motions, due 
to some processes (like magnetic reconnection) localized 
deep inside the magnetotail. On the basis of CLUSTER 
observations of reconnection events, a relationship be- 
tween the flapping motion and the reconnection process 
was investigated by 8]. During the reconnection events 
the current sheet exhibits strong flapping motions that 
propagate towards the flank of the tail. 

With regard to a theoretical aspect of the problem, 
the Ballooning-type mode in the curved current sheet 
magnetic field was claimed to be able to propagate az- 
imuthally in flankward directions from the source 
This ballooning theory was applied in the WKB approxi- 
mation implying the condition that the wave length scale 
is much less than the curvature radius. This condition 
can hardly be fulflUed in the plasma sheet with a small 
normal component of the magnetic field. Another point 
is that according to the theory of 0, both "kink" -like 
and "sausage" -like deformations of the current sheet are 
equally possible, and the question arises about a reason, 
why the observed flapping perturbations of the current 
sheet are mainly associated with the "kink" -like wave 
modes. 

In this paper, we propose a new approach to explain 
the existence of the "kink" -like flapping wave oscillations 
propagating across the current sheet. In a framework of 
a rather simple magnetohydrodynamic consideration, we 
elucidate a physical reason of the flapping wave oscilla- 
tions of the current sheet, which is related with gradients 
of the tangential and normal magnetic field components 
with respect to the normal and tangential directions, re- 
spectively. 
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STATEMENT OF PROBLEM 

A geometrical situation of the problem and coordinate 
system are illustrated in Fig. 1. We apply a system of 
incompressible ideal magnetohydrodynamics for nonsta- 
tionary variations of plasma sheet parameters 
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dt 
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ot 

V-V = 0, V-B = 0. 
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(2) 
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Here V, B,p, P are the velocity, magnetic field, density 
and total pressure, respectively. The total pressure is 
defined as the sum of the magnetic and plasma pres- 
sures. We consider specific wave perturbations prop- 
agating across the magnetic field lines, which are much 
slower than the magnetosonic modes. In this case the 
incompressible approximation seems to be appropriate. 

We focus our study on the very slow wave modes ex- 
isting only in the presence of a gradient of the com- 
ponent in the magnetotail current sheet along the x di- 
rection. The background conditions are considered to be 
rather simple with a slow dependence of the compo- 
nent on the X coordinate 



B,=B*b,iz), B, = eB*h,{x), By ^ 0, 
V = 0, y^y/A, z = z/A, x^x/L^. 



(4) 



Here A is a thickness of the current sheet, and is a 
length scale of the B^ variation along the current sheet. 

We introduce normalized small perturbations marked 
by sign "tilde" which are considered to be functions of 
time and two spatial coordinates (y, z) 

B^^B*iK + b^iz)), By^eB*by, 

B, = eB*ib,ix) + b,), P = Po + PP*V(4^), 

Vx = V^Va, Vy = VyVA, Vz = VyVA, 

i^tVA/A, Va=B*/^A^\ v^A/Lx. (5) 

Here Pq is the background total pressure, the parameter 
£ means the ratio of the background normal and max- 
imal tangential components of the magnetic field, and 
the parameter v characterizes the gradient of the nor- 
mal magnetic field component. For the background con- 
ditions considered in our model {Bz{x)^ Bx{z)), equation 
V-B — is fulfilled for arbitrary independent parameters 
e and u. 

Linearizing Eqs. ([TH3]) for the normalized perturba- 
tions, neglecting high order terms ~ v^e, and ^ e^, we 
assume v ^ e and retain the main term ^ ve. 

Substituting Fourier harmonics (oc exp{iuji — iky)) , we 
obtain finally a system of equations for Fourier ampli- 
tudes 
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T~ dP ~ dbz , ^ 

iujVy — ikP — 0, iujVz H — 7^ = svbx—;^, (7) 

dz dx 

-I , dvz ^ ^ dbz ._- dvy 

iLobz — bz—r:r + vVx—r:r = 0, lujby — bz—f^ — 0, (8) 
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0. (9) 



In this system of equations the derivative dbz/dx is as- 
sumed to be constant, and all other quantities are consid- 
ered to be not dependent on the x coordinate. Therefore 
Eqs. (l6]-l9]) are treated as a system of ordinary differential 
equations with respect to the z coordinate. Excluding bx 
and bz in Eq. ([6|), we derive 

/I , ^^ , ^ dbx db, 

Vx{~u;^ + U{z))=0, U{z)^ei^^^. (10) 

Generally, for a nonconstant U{z), Eq. (fTO|) yields Vx = 0. 

From Eqs. (l7]-[9]), we finally obtain a second order ordi- 
nary differential equation for the Vz velocity perturbation 

Further for simplicity we consider a piecewise constant 
function U (z) 

U{z) = ei^, ~1 <z <1; C/(z) = 0, |z| > 1, (12) 

which means that the current density is assumed to be 
constant within the current sheet. 



RESULTS 

A choice of the piecewise constant function U (z) allows 
us to find analytical solutions which are of two kinds, 
"kink" -like and "sausage" -like modes. The "kink" -like 
mode is characterized by displacement of the current 
sheet center, and even function Vz{z) 

Vz^Cexp{-k{\z\-l)), |z| > 1; (13) 
d^=Dcos{Xz), A = k.Jev/u}'^ ~ 1, |z| < 1. (14) 

An odd function Vz{z) is relevant to the "sausage" -like 
mode characterized by variations of the thickness of the 
current layer without a displacement of its center 

Vz = C exp(— fc(z — 1)), z > 1; 
Vz — ~C exp(k{z + 1)), z < —1; 
iz=Dsm{Xz), A = kyJevjuP- - 1, |z| < 1. (15) 

Applying continuity conditions for Vz and the first deriva- 
tive dvz/dz at the current layer boundaries, we obtain 
algebraic system corresponding to the "kink" mode 



C = i:'cos(A), fcC = Ai:'sin(A), 



(16) 
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and also we find a system for the "sausage" mode 
C = D sin(A), -kC = \D cos(A). 



(17) 



Setting the determinants to vanish, we derive two equa- 
tions corresponding to the "kink" and "sausage" modes, 
respectively 

tan(A) = ^ ("kink"); tan(A) = ("sausage"). (18) 
A k 

These equations have discrete sequences of roots 
Ai, A2, ...A„, ... . The main root is the minimal A which 
corresponds to the maximal frequency. 

By numerical solving these equations, we obtain two 
main roots Xk,s which determine the dimensional frequen- 
cies ujk,s as functions of wave number for the "kink" and 
"sausage" modes 



kA 



fc2A2 + XI , 



UJf = 



I 1 dB^dB^ 
Anp dz dx 



(19) 



Here ujf means a characteristic flapping frequency pro- 
portional to the square root of the multiplication of 
two gradients of the background magnetic field compo- 
nents, dBx/dz and dBz/dx. The dimensionless functions 
^k,s/^f are presented at the top panel in Fig. 2. Fre- 
quencies are monotonic functions of wave number, and 
they increase to the maximal asymptotic value ojf for 
fcA 00. The group wave velocity is shown in Fig. 2 
as functions of wave number (the second panel). It de- 
creases monotonically to zero for increasing wave num- 
bers. 

The flapping wave perturbations become unstable 
when the multiplication of two magnetic gradients be- 
comes negative. In particular, for the Earth's plasma 
sheet this condition corresponds to the case of decreas- 
ing Bz component towards Earth. The growth times of 
the instability for the "kink" and "sausage" modes are 
given by formulas 



Tk,s = Tf- 



fcA 



I -1 dB., dB, 

Anp dz dx 



.(20) 



The instability growth times {Tk,s/Tf) are shown in Fig. 2 
(bottom panel) as functions of wave number for the two 
wave modes. One can see from the figure that the unsta- 
ble "kink" mode develops much faster than the sausage 
mode. In particular, for kA = 0.7 the ratio of growth 
times is Ts/rk = 2. Fig. 3 illustrates a perturbation of 
the current sheet and the directions of plasma motion 
corresponding to the "kink" mode flapping. 

A qualitative explanation of the flapping instability 
and waves corresponding to the obtained solution is the 
following. Let us consider a plasma element of a unit vol- 
ume at the center of the current layer as shown in Fig. 4. 



Along the z direction the resulting force acting on this 
plasma element is a difference of two forces caused by the 
magnetic stress and the total pressure gradient. In equi- 
librium state, the resulting force Fz vanishes, and the 
total pressure gradient compensates the magnetic stress 



dP ^ dB^ 
dz Att ^ dx 



(21) 



In the new position of the magnetic tube element, the 
resulting force will be 

1 ,dBz 1 , fdB^dBz\ , , 

This force accelerates plasma in the z direction 



P- 



dHz_ 
W 



1 dB^ dBz 
4tt dz dx 



(23) 



This equation yields the characteristic flapping frequency 
Lu f which is proportional to the square root of the gradi- 
ents of the magnetic field components. This qualitative 
explanation of the instability is illustrated in Fig. 4 where 
panels (a) and (b) correspond to the stable and unstable 
situations, respectively. 

For example, we estimate this frequency for the pa- 
rameters which seem to be reasonable for the conditions 
of the current sheet in the Earth's magnetotail, 

B^ = 20 nT, Bz = 2nT, A - Re, Up = 0.1cm"^ 

fcA = 0.7, dBjdx ~ BJL,, ^ 5Re. (24) 

For these parameters we find the characteristic flapping 
frequency ujf ^ 0.03 s~^, and also the group velocity 
Vg = 60km/s. 



SUMMARY 

The flapping instability and waves are analyzed for a 
current sheet in a presence of two gradients of the B^ and 
Bz magnetic field components along the z and x direc- 
tions, respectively. These both gradients play a crucial 
role for the stability of the current sheet. The instability 
occurs in the regions of the current layer where the mul- 
tiplication of two gradients is negative. In particular, 
the instability can arise in a vicinity of a localized thin- 
ning of the current sheet (Fig. 4b). In stable regions, the 
flapping waves are associated with the so called "Bursty 
Bulk Flows" or BBF's which are the magnetic tubes 
rapidly moving through the center of the current sheet 
towards the Earth. These BBF's are considered to be 
the sources of the flapping wave oscillations propagating 
from the center of the current sheet towards the flanks 
in the ±y directions. 

The analytical solution is obtained for the simplified 
model of the current layer with a constant current den- 
sity. The frequency and the growth rate for the "kink" 
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FIG. 1: Geometrical situation of the problem 
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FIG. 2: Frequency, group velocity, and instability growth time 
as functions of wave number for two wave modes 



pronounced, because the flapping waves propagate in the 
direction perpendicular to the plane of the background 
magnetic field lines. The magnetic field planes are just 
shifting with respect to each other. 




® 



FIG. 3: Illustration to the "kink" mode. Perturbation of 
the current sheet and the corresponding directions of plasma 
motion 
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FIG. 4: Illustration to the "kink" flapping waves (a) and in- 
stability (b) in cases of positive and negative gradient of B^. 
Displacements of the magnetic tubes are shown. 



mode are found to be much larger than those for the 
"sausage" mode. For both modes, the frequencies are 
monotonic increasing functions of the wave number. The 
corresponding wave group velocities are decreasing func- 
tions of the wave number, and they vanishes asymptoti- 
cally for high wave numbers. 

For the typical parameters of the Earth's current sheet, 
the group velocity of the "kink" -like mode is estimated as 
a few tens of kilometers per second that is in good agree- 
ment with the CLUSTER observations. A strong de- 
crease of the group velocity for high wave numbers means 
that the small scale oscillations propagate much slower 
than the large scale oscillations. Because of that, the 
propagating flapping pulse is expected to have a smooth 
gradual front side part, and a small scale oscillating back- 
side part. 

The neglected second order terms 0(£^) are responsi- 
ble for the small effects related to the Alfven waves prop- 
agating in the z direction. These second order effects are 
subjects for future study. For the double gradient flap- 
ping waves studied in our model, magnetic tension is not 
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